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Abstract 

Inclusive Standard Model Higgs boson pair production at hadron colliders 
has the capability to determine the Higgs boson self-coupling, A. We present a 
detailed analysis of the gg —»• HH —► ( W + W~){W + W~ ) —> {jj£ ± v)(jjt' ± v) 
and gg —> HH —> (W+W~){W+W~) -> {jj£±v^vi'^v) (£, l', l" = e, g) 
signal channels, and the relevant background processes, for the CERN Large 
Hadron Collider, and a future Very Large Hadron Collider operating at a 
center-of-mass energy of 200 TeV. We also derive quantitative sensitivity lim¬ 
its for A. We find that it should be possible at the LHC with design luminosity 
to establish that the Standard Model Higgs boson has a non-zero self-coupling 
and that A/A sm can be restricted to a range of 0 - 3.8 at 95% confidence level 
(CL) if its mass is between 150 and 200 GeV. At a 200 TeV collider with an 
integrated luminosity of 300 fb -1 , A can be determined with an accuracy of 
8 — 25% at 95% CL in the same mass range. 
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I. INTRODUCTION 


The CERN Large Hadron Collider (LHC) is widely regarded as capable of directly observ¬ 
ing the agent responsible for electroweak symmetry breaking and fermion mass generation. 
This is generally believed to be a light Higgs boson with mass run < 200 GeV [1], a mass 
region for which discovery is covered by multiple channels, most notably by the decay that 
will be of interest here, H —> W + W~ [2], Once a Higgs boson candidate has been observed, 
emphasis will shift to a precise determination of its properties. The LHC promises complete 
coverage of Higgs decay scenarios [3], including general parameterizations in the Minimal 
Supersymmetric Standard Model (MSSM) [3,4], invisible Higgs decays [5], and possibly even 
Higgs boson decays to muons [6]. With mild theoretical assumptions and an integrated lu¬ 
minosity of 200 flu 1 , the Higgs boson total width, V H , and the gauge and various Yukawa 
couplings can be determined [7-9] with a precision of 10 — 30% [10]. At an e + e - linear col¬ 
lider with a center of mass energy of 350 GeV or more, these measurements can be improved 
by up to a factor 10 [11], if an integrated luminosity of 500 flu 1 can be achieved. 

While these studies have shown that future colliders promise broad and significant ca¬ 
pability to measure various properties of the Higgs sector, what remains is to determine the 
actual Higgs potential. This appears in the Lagrangian as 

V(4>) = _Au 2 ($ f $) + A(4> f 4>) 2 , (1) 

where is the Higgs field, v = (v^Gf) - 1 / 2 is the vacuum expectation value, and Gf is the 
Fermi constant. In the Standard Model (SM), 

A = A 5M = ||. (2) 

Regarding the SM as an effective theory, the Higgs boson self-coupling A is per se a free 
parameter. S'-matrix unitarity constrains A to A < 87t/3 [12], An anomalous Higgs boson 
self-coupling appears in various beyond the SM scenarios, such as models with a composite 
Higgs boson [13], or in two Higgs doublet models, for example the MSSM [14], To measure 
A, and thus determine the Higgs potential, at a minimum experiments must observe Higgs 
boson pair production. 

Strictly speaking, both the trilinear Higgs boson coupling gHHH and the quartic coupling 
Qhhhh have to measured separately in order to fully determine the Higgs potential. While 
9hhh can be measured in Higgs pair production, triple Higgs production is needed to probe 
Qhhhh- Since the cross sections for HHH production processes are more than a factor 10 3 
smaller than those for Higgs pair production at linear colliders [15], and about an order 
of magnitude smaller at hadron colliders [16], the quartic Higgs boson coupling will likely 
remain elusive even at the highest collider energies and luminosities considered so far. In 
the following we therefore restrict ourselves to gHHH which is related to A by gHHH = 3Au. 

Several studies of Higgs pair production in e + e~ collisions have been conducted over 
the past few years [15,17-19], and quantitative sensitivity limits for A have been derived 
for several proposed linear colliders with center of mass energies spanning the range from 

= 500 GeV to 3 TeV. For example, a study employing neural net techniques found that 
A can be measured with a precision of about 20% at a linear collider with a/s = 500 GeV 
and an integrated luminosity of 1 ab _1 , if m# = 120 GeV [18]. In contrast, the potential 
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of the LHC to probe the Higgs boson self coupling has begun to be explored only recently. 
A survey of Higgs pair production and background processes at an upgraded LHC, which 
would gather 20 times the amount of data expected in the first run (dubbed SLHC), was 
presented in Ref. [20]. In Ref. [21], we discussed the prospects for determining A at the LHC 
with design luminosity in Higgs pair production via gluon fusion and subsequent decay to 
same-sign dileptons via weak gauge W ± bosons, 

gg ^ hh ^ (w + w~)(w + w~) - (3) 

where £ and £' are any combination of electrons and muons, and presented quantitative 
estimates of sensitivity limits for A for 150 GeV < mn < 200 GeV and various integrated 
luminosities. Finally, in Ref. [22], a pythia [23] based study of several final states resulting 
from gg —> HH —> ( W + W~){W + W ~) and gg —> HH —> ( W + W~){ZZ ) for LHC and SLHC 
was conducted. 

In this paper we present a more detailed and extended analysis of Higgs pair produc¬ 
tion via gluon fusion at hadron colliders. In Ref. [21], we included only the two largest 
sources of background, W ± W + W~jj and ttW ± production, in our calculation. The effect 
of the remaining background contributions on the sensitivity limits for A was estimated 
by scaling the combined WWWjj and ttW cross section by a factor 1.1, as suggested by 
Ref. [20], which found that the remaining backgrounds are small. Here we present a more 
complete calculation of the background which includes, in addition to WWWjj and ttW 
production, W ± W ± jjjj, W ± Zjjjj , ttZ, ttj, tttt, W + W~W + W~ and WWZjj production. 
Furthermore, we discuss the potential size of backgrounds arising from overlapping events 
and double-parton scattering. We also extend our previous analysis by considering the three 
lepton channel 

gg —> HH —> ( W + W~)(W + W~) -> (4) 

(£, £', l" = e, /i), and the background processes which affect it. In addition to the LHC and 
SLHC, we calculate signal and background cross sections, and derive sensitivity bounds for A, 
for a Very Large Hadron Collider (VLHC), assuming a pp collider operating at y/s = 200 TeV 
with a luminosity of C = 2 x 10 34 cm _2 s _1 . These parameters correspond to one of the 
options listed in Ref. [24], 

The remainder of this paper is organized as follows. In Secs. II and III we outline our 
calculation of signal and background processes for the same-sign dilepton and three lepton 
final states, respectively. In Sec. IV, we derive sensitivity limits for A for various integrated 
luminosities at the LHC and VLHC. Our conclusions are given in Sec. V. 


II. THE SAME-SIGN DILEPTON FINAL STATE 

There are several mechanisms for pair production of Higgs bosons in hadronic collisions. 
Higgs boson pairs can be produced via gluon fusion, gg —> HH [16,25], weak boson fusion, 
qq —> qqHH [26], associated production with W or Z bosons, qq —> VHH, V = W, Z [27], 
and associated production with tt pairs, gg, qq —> ttHH [20]. At the LHC, inclusive Higgs 
boson pair production is dominated by gluon fusion. The weak boson fusion process, and 
associated production with W/Z bosons or tt pairs yield cross sections which are about a 
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FIG. 1. Representative Feynman diagrams for the process gg —> HH. 


factor 10 and 30 smaller than that for gg —» HH [25,20]. Since Higgs pair production at the 
LHC is rate limited, we concentrate on the gluon fusion process in the following. 

For rtiH < 140 GeV, the dominant decay mode of the SM Higgs boson is H —> bb, and 
the QCD bbbb background overwhelms the gg —> HH signal [28]. For m# > 140 GeV, 
H —> W + W~ dominates, and the W + W~W + W~ final state has the largest individual 
branching ratio. If all W bosons decay hadronically, QCD multi-jet production dwarfs 
the signal. A similar result is obtained for the £±v + 6 jet (only one W boson decays 
leptonically), and + 4 jet (one W + W~ pair decays leptonically) final states, where 

W+ multi-jet and W + W~+ multi-jet production provide very large backgrounds. This 
leaves the same-sign dilepton final states, {jj£ ± v){jj£' ± v), modes where three W bosons 
decay leptonically and one decays hadronically, and the all-leptonic decay modes. The 
latter suffer from a large suppression due to the small WWWW —> M + 4z/ branching ratio 
of (0.216) 4 = 0.0022 (BR(kF —> Iv) = 0.216, £ — e, ji). In the following we therefore only 
consider the {jj£ ± v){jj£ t± v) and (jj£ ± v)(£ ,± v£" T v) hnal states. 

In this section we discuss in detail the calculation of signal and background cross sections 
for the (jj£ ± ^')(jj£ ,± ^) hnal state. The three lepton hnal state will be considered in Sec. III. 


A. Calculation of the signal cross section 

The Feynman diagrams contributing to gg —» HH in the SM consist of fermion triangle 
and box diagrams (see Fig. 1) [16]. Non-standard Higgs boson self-couplings only affect 
the triangle diagrams with a Higgs boson exchanged in the s-channel. We calculate the 
gg —> HH —> (W + W~)(W + W~) —> {jj£ ± v)(jj£' ± v) cross section using exact loop matrix 
elements [16]. As demonstrated in Ref. [21], the infinite top quark mass limit, which is 
commonly used in place of exact matrix elements to speed up the calculation, reproduces 
the correct total cross section for HH production to within 10% to 30% for Higgs masses 
between 140 GeV and 200 GeV, but produces completely incorrect kinematic distributions. 
The intermediate Higgs and W bosons are treated off-shell using finite widths in the double 
pole approximation in our calculation. Decay correlations for the H —> W + W~ —> 4 fermion 
decays are fully taken into account [29]. 

Signal results are computed consistently to leading order QCD with the top quark mass 
set to m t = 175 GeV and SM HWW and top quark Yukawa couplings, and the renormaliza¬ 
tion and factorization scales are taken to be the Higgs boson mass [16]. The contributions 
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of all other quarks to the gg —» HH box and triangular diagrams are strongly suppressed 
due to their Hff couplings, which are proportional to the fermion mass. The effects of 
next-to-leading (NLO) order QCD corrections are included in our calculation by multiply¬ 
ing the differential cross section by an overall factor K = 1.65 at LHC and K = 1.35 at 
VLHC energies (K- -factor) for scale choice g = mn, as suggested by Ref. [30] where the 
QCD corrections for gg — > HH have been computed in the large rn t limit. Although this 
approximation cannot replace a calculation of the full NLO QCD corrections to gg —> HH, 
it is expected to work well in this particular case: it is well known from single Higgs bo¬ 
son production via gluon fusion [31] that the two-loop QCD corrections to the one-loop 
Higgs production amplitude are well approximated by multiplying the leading order one- 
loop cross section for a finite top quark mass by the AT-factor obtained in the large rnt limit. 
This feature can be easily understood by recalling that the dominant corrections originate 
from radiation off the initial state gluons, which is universal. 

The lowest order gg —> HH cross section exhibits a rather strong dependence on the 
renormalization and factorization scales. For example, for rrifj = 200 GeV and g 2 = rn 2 H , 
and using exact loop matrix elements, one obtains a total cross section of 8.26 fb. For 
g 2 = s, on the other hand, one finds a cross section which is almost a factor 1.5 smaller. 
NLO QCD corrections, albeit computed in the large m t limit, significantly reduce the g 
dependence: varying the scale g from g 2 = m 2 H to g 2 = s, the cross section including the 
A'-factor decreases by a factor 1.25 instead of 1.5 without the A'-factor. 

For g 2 = rnjj , the A'-factor at VLHC energies is smaller than that obtained for the 
LHC. At higher energies, smaller parton momentum fractions are probed. This results in an 
increased sensitivity of the cross section to the choice of factorization scale which partially 
compensates the variation of the cross section with the renormalization scale. 

In all our calculations we use a value for the strong coupling constant of a s (M z ) = 
0.1185 [32], All signal and background cross sections are computed using CTEQ4L [33] 
parton distribution functions. 

The kinematic acceptance cuts for both signal and backgrounds in the (jjt'D/) 
channel at the LHC and VLHC are: 

p T (j) > 30, 30, 20, 20 GeV, p T (i) > 15, 15 GeV, (5) 

\v(j)\ < 3-0, \v(£)\ < 2.5, (6) 

AR(jj) > 0.6, A R{j£) > 0.4, A R(££) > 0.2, (7) 

2 21 1 /^ 

where A R = (A0)" + (A//)N is the separation in the pseudorapidity - azimuthal angle 
plane. In addition we require the four jets to combine into two pseudo-IF pairs with invariant 
masses 

50 GeV < m(jj) < 110 GeV, (8) 

and assume that this captures 100% of the signal and backgrounds. We do not impose a 
missing transverse momentum cut which would remove a considerable fraction of the signal 
events. Detector resolution effects are not taken into account in our calculation. 

As we have shown in Ref. [21], at the LHC, the main background from WWWjj and 
ttW production can be reduced by about 45%, with little impact on the signal, by imposing 
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FIG. 2. The pp —> HH —► (IV + IV _ )(fF + IV _ ) —> + 4j differential cross section as a 

function of the minimum lepton transverse momentum, PTmin{£)i for pp collisions at y/s = 14 TeV. 
Results are shown for rrifj = 150 GeV (solid line), mn = 160 GeV (dashed line), m# = 180 GeV 
(dotted line), and mu = 200 GeV (dash-dotted line). 


a more restrictive jet-jet separation cut of A R(jj) > 1.0. In contrast, at VLHC energies, 
there is little gain in tightening the A R(jj) cut. Subsequently, we therefore require 

A R(jj) > 1.0 at the LHC, and (9) 

A R{jj) > 0.6 at the VLHC (10) 

in the (jjWz) channel. 

Our choice of p T cuts for jets and leptons in Eqs. (5) - (7) is driven by the goal of retaining 
as much signal as possible while ensuring that the LHC experiments, ATLAS and CMS, can 
record + 4 j events when operating at the LHC design luminosity of £ = 10 34 cm -2 s _1 . 

Figures 2 and 3 show the pp —> HH —> {W + W~){W + W~) —> £ ± £ /± + 4 j differential cross 
section at the LHC as a function of the lepton and jet minimum transverse momentum, 
respectively. Qualitatively similar results are obtained at VLHC energies. The differential 
cross section peaks at low values of PTmin and falls very quickly with increasing values of 
the minimum transverse momentum, in particular in the jet case. In order to maximize the 
signal cross section, the lepton and jet Pt thresholds are thus chosen as low as possible and 
yet be compatible with the requirements of ATLAS and CMS to successfully record such 
events. 

We shall also use the cuts listed in Eqs. (5) - (7) for a luminosity upgraded LHC operating 
at £ = 10 35 cm -2 s _1 , and at the VLHC. Preliminary studies have concluded [20,34] that cuts 
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FIG. 3. The pp —► HH —> ( W + W~)(W + W ~) —> ^ ± ^ /± + 4j differential cross section as a 
function of the minimum jet transverse momentum, PTmmij): for pp collisions at yfs = 14 TeV. 
Results are shown for ran = 150 GeV (solid line), mjj = 160 GeV (dashed line), m# = 180 GeV 
(dotted line), and mu = 200 GeV (dash-dotted line). 


similar to those listed in Eq. (5) - (7) should be sufficient, although increased background 
from event pilcup is expected to degrade detector performance, in particular at the SLHC. 


B. Calculation of the backgrounds 

The SM backgrounds of interest are those that produce two same-sign leptons and four 
well-separated jets which reconstruct in two pairs to a window around the W boson mass. 
The largest contribution originates from W ± W + W~jj production, followed by ttW ± where 
one top quark decays lcptonically, the other hadronically, and neither b quark jet is tagged. 
Other backgrounds which contribute are: W ± W ± jjjj production; tttt production, where 
none of the b quark jets are tagged, and additional jets or leptons are not observed; W ± Zjjjj, 
tiZ and W + W~ Zjj production with leptonic Z decay (including off-shell photon interfer¬ 
ence) where one lepton is not observed; and ttj events where one b quark decays semilep- 
tonically with good hadronic isolation, and the other is not tagged. In addition, in a high 
luminosity environment, one has to worry about backgrounds from overlapping events and 
double parton scattering. 6-quarks are assumed to be tagged with an efficiency of 50% 
throughout. We do not apply an explicit A'-factor for the backgrounds here; however, we do 
later include the potential effect of QCD corrections on the backgrounds when we extract 
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limits on the Higgs self-coupling (see Sec. IV). 


1. The WWWjj, ttW, W + W W + W and tttt backgrounds 

We simulate these backgrounds at the parton level using exact matrix elements generated 
with madgraph [35]. The WWW jj background has a significant contribution from WH{-^ 
W + W~)jj production. For WWW jj production we evaluate the strong coupling constant 
a s and the parton distribution functions at a scale p given by /a 2 = YPtj where the sum 
extends over all final state particles; for ttW production we take p = m t + M w / 2, and in the 
tttt case we use p = 2 m t . Top quarks are generated on shell (narrow width approximation), 
while all W bosons are allowed to be off shell. Events with one or more tagged h quarks are 
rejected. In the tttt case we merge jets if their separation in the pseudorapidity - azimuthal 
angle plane is AR(jj) < 0.6. The W + W~W + W~ background has a significant contribution 
from WWH (-> WW) production. The W + W~W + W~ cross section at the LHC (VLHC) is 
found to be a factor 5 to 25 (25 to 80) smaller than that of WWW jj production, depending 
on the mass of the Higgs boson. In this analysis we therefore ignore the W + W~W + W~ 
background. 


2. The W^W^jjjj background 

Although MADGRAPH is able to generate exact matrix elements for W ± W ± jjjj pro¬ 
duction, the large number of contributing Feynman diagrams (more than 6,000) makes a 
full matrix element based calculation of the W ± W ± jjjj background impractical. In order 
to estimate the W ± W ± jjjj cross section we thus have interfaced the matrix elements for 
pp —> W ± W ± jj [36] with pythia, which produces the two additional jets in a leading-log 
shower approximation. The strong coupling constant a s and the parton distribution func¬ 
tions are evaluated at a scale p given by p 2 = YPti where the sum extends over all final 
state particles. 


3. The ttZ, W ± Zjjjj and WWZjj backgrounds 

We calculate the ttZ cross section using exact matrix elements. Since more than 15,000 
Feynman diagrams contribute to W ± Zjjjj production, we estimate the cross section for this 
process by interfacing the W ± Zjj matrix elements [37] with PYTHIA, similar to W ± W ± jjjj 
production. Off-shell photon interference effects are taken into account in both cases. Both 
processes contribute to the background only if one of the leptons from Zj off-shell photon 
decay is missed. We consider a lepton to be missed if it has pt < 10 GeV or \p\ > 2.5. If the 
lepton is within a cone of A R < 0.2 from a detected lepton and has 1 GeV < p T < 10 GeV, 
then the detected lepton is not considered isolated and the event is rejected. The strong 
coupling constant and the parton distribution functions are evaluated in ttZ ( W ± Zjjjj ) 
production at p — m t + Mz/2 (p 2 = YPt )• I n order to avoid the collinear singularity 
when the missed lepton is collinear with an observed lepton (which is only relevant if the 
missed lepton has Pt < 1 GeV), finite lepton masses must be maintained in the calculation 
of the W ± Zjjjj and ttZ processes [9], the latter also including additional considerations to 



maintain gauge invariance of the calculation, clue to finite top quark width which must be 
included for these events. ttZ events are rejected with a factor 4, which approximates the 
fraction of events with one or more tagged 6 quarks. 

The size of the WWZjj background can be estimated from the ratio of the WWZ and 
WWW cross sections, together with the suppression factor which arises from requiring that 
one lepton is missed, and the WWW jj rate. We find that the WW Zjj cross section is 
about a factor 30 to 60 smaller than the WWW jj cross section, depending on the mass of 
the Higgs boson. We thus ignore the WW Zjj background in the following. 


4 • The ttj background 

We calculate ttj —> W + bW~bj production where one of the 6-quarks decays semilcptoni- 
cally with an isolated lepton using exact tt+parton and b —» civ matrix elements. The lepton 
from 6-decay is considered not isolated if the charm quark is within a cone of A R < 0.4 
from the lepton and has px(c) > 3 GeV. Events are rejected with a factor two, which ap¬ 
proximates the fraction where the second 6 quark would be tagged. In order to regularize 
the soft parton px distribution so as to reproduce the px distribution of the hard recoil 
system (tt) from a full resummation calculation [38], while preserving the normalization 
of the hard tt cross section, we use the truncated shower approximation (TSA) [39]. The 
advantage is that QCD matrix elements at tree level contain the full information on angular 
distributions and hardness of additional jet emission. A parton shower approach as used 
in Ref. [22] would not immediately give reliable answers unless both color coherence and 
the choice of scale are implemented correctly, matching the answer given by QCD matrix 
elements for sufficiently hard partons. In practice this is achieved by integrating the full 
tree-level tt+parton matrix elements over phase space down to px (parton) > 1 GeV, and 
multiplying the result by a factor 1 — exp[— pf.(parton) /pxsAi w h ere Ptsa is adjusted to 
achieve a t tj = a t f. For a scale choice of p 2 = rn 2 , we find Pxsa = 15, 50 GeV for the LHC 
and VLHC, respectively. This should not be construed as an attempt to mimic a full NLO 
calculation of ttj production, which is not available, but it does represent an improvement 
over using only PYTHIA tt matrix elements with additional partonic emission from shower¬ 
ing, and includes a well-motivated approach to controlling the soft singularity present in the 
exact O(a^) matrix elements. 


5. Other backgrounds: overlapping events and double parton scattering 

At a high-luminosity intersection point of colliding beams, more than one event may 
occur per bunch crossing. In principle, such overlapping events can be recognized by a total 
visible energy measurement or by tracing some final particle tracks back to distinct event 
vertices, but in practice this may not always be possible. In this case, the overlap of, e.g, two 
single H —> W + W~ events, a W + W~ and a Wjj event, or a tt and a single W event may 
constitute a potentially dangerous source of background for the Higgs boson pair production 
signal. If the vertices of the overlapping events cannot be resolved, the effective cross section 
for overlapping events is given by [36] 
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TABLE I. Luminosity, bunch spacing, and luminosity per bunch crossing for the LHC [20], 
SLHC [20] and VLHC [24], 


machine 

luminosity C 

bunch spacing Ar 

Ac 

LHC 

10 34 cm" 2 s" 1 

25 ns 

7^ 

o 

4 

SLHC 

10 35 cm" 2 s" 1 

12.5 ns 

7 

S 

GO 

O 

VLHC 

2 x 10 34 cm" 2 s” 1 

18.8 ns 

To 

4 


<Jov(l,2) = ^<7(1)<7(2) Ac, (11) 

where cr(l) and cr(2) are the cross sections for the two overlapping processes, and C bc is the 
luminosity per bunch crossing. It is given by 

jC hc = jCAr, ( 12 ) 

where C is the instantaneous luminosity and At is the bunch spacing. The values for C, At 
and C bc at the LHC, SLHC and VLHC are listed in Table I. 

The P ± £' ± jjjj final state can also be produced via the independent scattering of two 
pairs of partons in the incident protons. The cross section for double parton scattering is 
given by Eq. (11) with the factor C bc /2 replaced by l/cr e ff. The parameter a e s & 15 mb [40], 
the effective cross section, contains all the information about the non-perturbative structure 
of the proton in this simplified approach. It is believed that <r e ff is largely independent of 
the center of mass energy [41]. Comparison of a e g and the values for C bc listed in Table I 
shows that the double parton cross section is about a factor 2 to 10 smaller than that from 
overlapping events. 


C. Numerical results 

The total cross sections within cuts (see Eqs. (5) - (7)) for signal and background pro¬ 
cesses at the LHC and VLHC are listed in Table II. At the LHC, with 300 fb ', at most 
about 50 signal events are produced. Outside of the Higgs boson mass range considered 
here, the number of signal events is too small to be useful. For m# < 150 GeV, this is due 
to the small H W*W branching ratio. For rrifj > 200 GeV, the gg —> HH cross section 
is too small. WWWjj and ttW production are the largest contributions to the background. 
The background from ttZ production where one of the leptons is lost is moderate. Although 
the cross section for WZjjjj production is substantial, this background can be separated 
rather easily from the signal as discussed below. WWjjjj and tttt production contribute 
negligibly to the background at the LHC. The tttt cross section is suppressed by the large 
top quark mass. The WWjjjj cross section is small because quark-gluon and gluon-gluon 
fusion processes do not contribute to same sign W pair production. 

The ttj cross section is extremely sensitive to the lepton pt cut imposed. Requiring 
Pt(£) > 15 GeV, the ttj background is of the same size as the ttZ background, and about 
a factor 3 smaller than the ttW background. Decreasing the p T (£) cut to 10 GeV, the ttj 
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TABLE II. Higgs pair signal and background cross sections (fb) for pp —» £^£ ,= + 4 j 
(£, i' = e, fi) at (a) the LHC (y% = 14 TeV) and (b) at the VLHC (y/s = 200 TeV), impos¬ 
ing the cuts listed in Eqs. (5) - (7), and as a function of the Higgs boson mass (GeV). The 
background labeled “pileup” represents a rough estimate of the combined WWWjj, tiW , ttZ, 
WZjjjj , WWjjjj and tttt cross section from overlapping events and double parton scattering. 
Cross sections at the SLHC are identical to those in the LHC case with the exception of the pileup 
cross section, which is about a factor 3.7 larger than at the LHC. The last column, labeled Btot , 
shows the total background cross section. 


m H 

HH 

WWWjj 

tiW 

ttZ 

(a) LHC 
tij 

WZjjjj 

WWjjjj 

tttt 

pileup 

Btot 


0.07 

0.36 

0.22 

0.05 

0.08 

0.15 

0.005 

0.002 

~ 0.03 

0.90 


0.19 

0.49 

0.22 

0.05 

0.08 

0.15 

0.005 

0.002 

~ 0.03 

1.03 


0.18 

0.40 

0.22 

0.05 

0.08 

0.15 

0.005 

0.002 

~ 0.03 

0.94 


0.08 

0.29 

0.22 

0.05 

0.08 

0.15 

0.005 

0.002 

~ 0.03 

0.83 






(b) VLHC 






m H 

HH 

WWWjj 

tiW 

tiZ 

tij 

WZjjjj 

WWjjjj 

tttt 

pileup 

Btot 

140 

2.2 

14.9 

5.8 

7.4 

7.7 

8.1 

0.13 

6.13 

~ 20 

70.2 

150 

6.5 

17.0 

5.8 

7.4 

7.7 

8.1 

0.13 

6.13 

~ 20 

72.3 

160 

15.8 

20.4 

5.8 

7.4 

7.7 

8.1 

0.13 

6.13 

~ 20 

75.7 

180 

16.0 

17.9 

5.8 

7.4 

7.7 

8.1 

0.13 

6.13 

~ 20 

73.2 

200 

8.7 

14.3 

5.8 

7.4 

7.7 

8.1 

0.13 

6.13 

~ 20 

69.6 


cross section increases by about a factor 10, overwhelming the Higgs pair signal 1 . On the 
other hand, if the minimum lepton transverse momentum is increased to 20 GeV which 
reduces the signal cross section by about 20%, the ttj background decreases by one order 
of magnitude and essentially becomes negligible. However, we emphasize that our matrix 
element based calculation of the ttj background should be viewed with some caution. Effects 
from hadronization, event pileup and extra jets from initial or final state radiation, as well as 
detector resolution effects may significantly affect the cross section. For a reliable estimate 
of the background, a full detector simulation, which is beyond the scope of this paper, is 
required. 

The lepton isolation requirement, together with the lepton px cut, the 6 —> civ branching 
ratio and the di-jet invariant mass cut suppress the ttj cross section by about a factor 10 6 . A 
similar suppression factor is also expected in £ ± ubb + ?jj production which also contributes to 
the background if one of the 6-quarks decays semileptonically and if the lepton from 6-decay 
is isolated. Using the result of Ref. [42], we estimate that the £ ± z/66 + 3 j cross section at the 
LHC is of 0( ICC 3 fb) which can safely be neglected. 


1 In Ref. [21] we required one lepton with a px > 10 GeV, and one lepton with px > 15 GeV. In 
this case, the ttj background is significantly larger than the signal (see also Ref. [22]). Increasing 
the lepton transverse momentum cut to 15 GeV solves this problem, while reducing the signal cross 
section by less than 10%. 
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Our numerical results for the overall normalization of the signal, the WWW jj, ttW , 
and the titi background processes agree reasonably well with those reported in Refs. [20] 
and [22]. For WZjjjj production, we find a cross section which is about a factor 10 larger. 
The discrepancy can be traced to the contribution from virtual photon exchange, which was 
not taken into account in [20,22], No result for ttZ production is given in Refs. [20] and [22], 
A meaningful comparison of our matrix element based calculation of the tij background and 
the pythia based estimate in [20,22] is not possible due to the strong dependence of the 
cross section on the lepton p T . 

Overlapping events and double parton scattering are not expected to contribute signifi¬ 
cantly to the background at the LHC. Contributions from these sources are listed in Table If 
in the column labeled as “pileup”. The numerical values listed were obtained by adding all 
overlapping event and double parton scattering contributions to the sources of background 
discussed in Sec. 11B, using Eq. (11) and the values of Cb c given in Table I. Since Eq. (11) 
assumes that that the vertices of the two overlapping events are not resolved, these values are 
likely conservative. For the SLHC, the pileup cross section in Table Ila has to be multiplied 
by a factor 3.7. Our results for the cross section from overlapping events and double parton 
scattering should be regarded only as order of magnitude estimates. Realistic simulations 
are needed to draw firm conclusions for this background. 

At a pp collider with a/s = 200 TeV, the cross sections of processes which are dominated 
by gluon fusion, such as the the gg —> HH signal, tttt, ttZ and tij production, are about 
a factor 100 - 3000 larger than that at the LHC. Since the cross section for titi production 
at the LHC is suppressed by the large invariant mass of four top quarks, the increase is 
particularly large for this process. In contrast, the cross sections of processes dominated 
by quark-gluon fusion or quark-quark scattering, such as WWWjj, tiW and WWjjjj 
production, increase by only a factor 25 - 45. As a result, the ttZ, tij and titi backgrounds 
are relatively more important at the VLHC. The cross sections due to overlapping events and 
double parton scattering increase by almost three orders of magnitude, and thus may well 
compete in size with WWW jj production, unless the vertex positions of the overlapping 
events are resolved. Since the signal is purely gluon induced, the overall signal to background 
ratio at the VLHC is about a factor 2 better than at the LHC. 

All the backgrounds are multi-body production processes, therefore the distribution of 
the invariant mass, V§, of the system peaks at values significantly above threshold. In 
contrast, the signal is a two-body production process for which the \/§ distribution will 
exhibit a sharper threshold behavior. Unfortunately, with two neutrinos present in the final 
state, \i§ cannot be reconstructed. However, we anticipate that the invariant mass of all 
observed final state leptons and jets given by (E) and p, are the energies and momenta of 
the jets and leptons) 

1 2 r -i 2 

m lis = J2 E i ~ J2 Pi ( 13 ) 

i=£,£', jets , jets 

will retain most of the expected behavior of the different production processes. Figures 4 
(LHC) and 5 (VLHC) clearly demonstrate that this is the case: the signal peaks at smaller 
values of m V i S than the background processes, especially for lower Higgs boson masses. 
This distribution, which was not considered in Ref. [20], is what makes possible a y 2 based 
test to improve extraction of the Higgs boson self-coupling (see Sec. IV). While detector 
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FIG. 4. Distribution of the invariant mass of the observable final state particles, m v i S , after all 
cuts, in pp —► ^ ± G ± + 4 j for the signal with a) rri jj = 150 GeV and b) run = 180 GeV, and all 
backgrounds (except for the contributions from overlapping events and double parton scattering) 
at the LHC. The dot-dashed curve shows the combined cross section of WZjjjj , WWjjjj and 
tttt production. 
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m V is (GeV) 

FIG. 5. Distribution of the invariant mass of the observable final state particles, m V i S , after all 
cuts, in pp —► ^ ± £ /± + 4 j for the signal (solid line) with a) mu = 150 GeV and b) run = 180 GeV, 
and all backgrounds (except for the contributions from overlapping events and double parton 
scattering) at the VLHC (dashed: WWWjj, dotted: ttW, long-dashed: ttZ, long-dash-dot: ttj). 
The dot-dashed curve shows the combined cross section of WZjjjj, WWjjjj and tttt production. 
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effects may smear out the tails of this distribution, or shift a peak slightly, it is a genuine, 
simple physics effect and care must be taken in any approximations used to simulate the 
backgrounds that this physics feature is retained. 

Since the WWWjj background has a significant contribution from WH{-^ W + W~)jj 
production, its rn ms distribution is similar in shape to that of the HH signal. As expected, 
the rn ms distributions of ttW, ttZ and ttj production peak at similar values, and are similar 
in shape. The dot-dashed lines in Figs. 4 and 5 represent the combined differential cross 
section of WZjjjj, WWjjjj and tttt production. It peaks at a much higher visible invariant 
mass than those of the other background processes. Whereas the signal is concentrated in the 
region m vis < 500 GeV, the background processes have a significant tail extending to m vis = 
1 TeV and beyond. This makes it possible to normalize the background using data from the 
m V i S > 500 GeV region. The simple procedure described in Sec. IIB 5 for estimating the cross 
section for overlapping events and double parton scattering is not suitable for calculating 
distributions. The m V i S distribution of these backgrounds therefore is not included in Figs. 4 
and 5. With increasing Higgs boson mass, the signal peak gradually moves to higher values of 
rn ms . The efficiency of the rn vls distribution as a discriminator thus decreases somewhat for 
niH > 180 GeV. Comparing the visible invariant mass distributions at the LHC and VLHC, 
the improved signal to background ratio at the VLHC in the region m vis = 200 — 500 GeV, 
where the Higgs pair signal is concentrated, is obvious. 

As noted before, all our calculations are consistently performed at leading order, ie. 
there are precisely four jets (partons) in the final state. In practice, one expects a significant 
fraction of the £ ± £ ,± + 4 j signal events to contain one or more extra jets originating from 
initial state gluon radiation. In such events, it is natural to construct rn ms from the four 
highest pt jets in the event. However, there is no guarantee that the extra jets are always 
the softest jets in the event. Since the m(jj ) requirement is rather loose (see Eq. (8)), it is 
conceivable that events where one (or more) of the four jets incorporated in rn ms originate 
from QCD bremsstrahlung. Hard QCD corrections could also lead to + 4j events 

where one of the jets from W decay does not pass the minimum pt cut for jets, but the 
additional bremsstrahlung jet does. Some of these events might also pass the m(jj ) cut. 
QCD corrections thus could affect the shape of the m V i S distribution. In order to draw firm 
conclusions, a full calculation of the NLO QCD corrections to gg —> HH with finite top 
quark mass effects is needed. Insight may also be gained from performing a calculation 
where the gg —> HH matrix elements are interfaced [43] with an event generator such as 
PYTHIA. 

In using pythia for the additional jet radiation, however, one has to be careful. As 
described previously, the radiation of soft and collinear jets from the initial state is the main 
source of the large (and top mass independent) QCD corrections to the total signal cross 
section. The initial state radiation modeled by pythia effectively resums the leading effects 
of precisely this radiation and includes it in the topology of the final state. Normalizing the 
rate to the leading order total cross section is therefore inconsistent and the result arbitrary 
(and not, as often is claimed, a conservative estimate), because the final state topology and 
the rate are computed in different approximations with a difference which is by no means a 
reduced higher order uncertainty. 

The effect of hard QCD corrections on the m V i S distribution may be reduced by limiting 
the number of possible 4 j combinations which satisfy the cut of Eq. (8). Approximately 
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60 — 65% (35 — 40%) of all signal events have one (two) 4 j combination satisfying Eq. (8); 
almost none have three 4 j combinations in the correct invariant mass range. When additional 
jets are present, many more combinations are possible. Adding the requirement that at most 
two 4 j combinations satisfy Eq. (8) may thus reduce the effect of hard QCD corrections on 
the m V i S distribution. It may also reduce the signal and background cross sections somewhat. 
The fraction of events where one or several QCD bremsstrahlung jets pass the cuts may also 
be reduced by shrinking the m(jj ) range in Eq. (8) (see also Ref. [22]). Our choice has been 
deliberately conservative. Reducing the di-jet invariant mass range to M\y ± 20 GeV may 
well be possible [44], This would also improve the signal to background ratio. 


III. THE THREE LEPTON FINAL STATE 

The calculation of signal and background cross sections for the (jj£ iz v)(£ ,± v£"^v) final 
state is similar to that described in Sec. II for the same sign di-lepton final state. Due to the 
smaller branching ratio for leptonic W decays, the cross section is expected to be somewhat 
smaller than that for the (jj£ ± v)(jj£ l± i') channel. The kinematic acceptance cuts for both 
signal and backgrounds in the (jj£ iz v)(£ ,± v£ ,/T v) final state are: 

pr(J) > 30, 20 GeV, Pt(£) > 15, 15, 15 GeV, (14) 

\VU )I <3-0, |/ 7 (0I < 2.5, (15) 

A R(jj) > 0.6, A R(j£) > 0.4, A R(U) > 0.2. (16) 

In addition, we impose the di-jet invariant mass cut of Eq. (8). 

Except for W ± W ± jjjj and W ± Zjjjj production, all processes discussed in Sec. IIB 
contribute to the background. For final states containing a same flavor opposite sign lep¬ 
ton pair, W ± £ + £~jj production constitutes an additional source of background. We have 
calculated the W ± £ + £~jj background using the exact matrix elements of Ref. [37]. The 
cross section for W ± £ + £~jj production is about three orders of magnitude larger than the 
Higgs boson pair signal, if photon exchange diagrams are taken into account. To increase 
the signal to background ratio, one can either impose a minimum £ + £~ invariant mass cut, 
or increase the A R(£ + £~) cut. Unfortunately, due to correlations between the momenta of 
the fermions in the H —» W + W~ — > 4 fermion decay [29], the £ + £~ invariant mass tends to 
be rather small in Higgs pair events. As a result, the signal to background ratio cannot be 
improved to better than 1:100 without reducing the signal cross section to an unacceptably 
low level. 

In the following we therefore only consider + 2 j, i % U, production. Excluding 

all final states containing a same flavor opposite sign lepton pair reduces the signal cross 
section by a factor 4. As a result, the cross section at the LHC becomes too small to be of 
interest; for an integrated luminosity of 300 fb -1 only 8 events are expected. We therefore 
present numerical results at VLHC energies only in this section. 

The total cross sections within cuts (see Eqs. (14) - (16)) for signal and background 
processes at the VLHC are listed in Table III. The signal cross section is about a factor 5 
smaller than in the same-sign dilepton case. The loss in signal however is at least partially 
compensated by the significantly improved signal to background ratio. This becomes more 
evident in Fig. 6 which shows the rri ms distribution for m# = 150 GeV and m# = 180 GeV. 
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TABLE III. Higgs pair signal and background cross sections (fb) for pp —> + 2 j 

(£, l' = e, p, t ^ £') at the VLHC (y% = 200 TeV), imposing the cuts listed in Eqs. (14) - (16), 
and as a function of the Higgs boson mass (GeV). The background labeled “pileup” represents a 
rough estimate of the combined WWWjj, ttW, ttZ , and tttt cross section from overlapping events 
and double parton scattering. The last column, labeled Btot, shows the total background cross 
section. 



HH 

WWWjj 

ttW 

ttZ 

tt.j 

tttt 

pileup 

Btot 


1.40 

1.43 

0.39 

1.39 

0.45 

0.47 

~ 2.4 

6.53 


3.06 

1.96 

0.39 

1.39 

0.45 

0.47 

~ 2.4 

7.06 


3.04 

1.71 

0.39 

1.39 

0.45 

0.47 

~ 2.4 

6.81 


1.66 

1.47 

0.39 

1.39 

0.45 

0.47 

~ 2.4 

6.57 


Here, rn ms is defined by Eq. (13), using the four-momentum vectors of the three charged 
leptons and the two jets. The largest contribution to the background in the region where the 
signal peaks comes from ttZ production. While the tij and tttt total cross sections are similar 
(see Table III), the tttt contribution to the background in this region is negligible. Most tttt 
events have visible invariant masses well in excess of 1 TeV. The signal to background ratio 
for the £ ± £ /=F £ /=F + 2 j channel is approximately a factor two better than for the same-sign 
dilepton final state (see Fig. 5). 

Since there are fewer jets present, initial state gluon radiation should have a smaller 
effect on the rn ms distribution in the three lepton final state than in the same-sign dilepton 
case. 


IV. DETERMINING THE HIGGS BOSON SELF-COUPLING 

The Feynman diagrams contributing to gg —» HH in the SM consist of fermion triangle 
and box diagrams [16]. Non-standard Higgs boson self-couplings affect only the triangle 
diagrams with a Higgs boson exchanged in the s-charmel. They contribute only to the 
J = 0 partial wave, and thus impact the rn ms distribution mostly at small values. This 
is illustrated in Fig. 7 for the t ± t' ± + 4 j final state with rrin = 180 GeV and two non¬ 
standard values of A hhh = A/A sm • Since box and triangle diagrams interfere destructively, 
the gg —> HH cross section for 1 < A hhh < 2.7 is smaller than in the SM. The absence of a 
Higgs boson self-coupling (A hhh = 0) results in a Higgs pair production cross section which 
is about a factor 3 larger than the SM result. Figure 7 also demonstrates that the m V i S 
distribution of the signal peaks at a smaller value than that of the combined background. 
This remains true for other Higgs boson masses, as long as rrin < 200 GeV. 

The shape change of the m vis distribution induced by non-standard values of A hhh can 
be used to derive quantitative sensitivity bounds on the Higgs boson self-coupling. We 
accomplish this by calculating 95% confidence level (CL) limits performing a y 2 test. The 
statistical significance is calculated by splitting the rn ms distribution into a number of bins, 
each with typically more than five events. In each bin the Poisson statistics are approximated 
by a Gaussian distribution. We impose the cuts described in Secs. II and III and combine 
channels with electrons and muons in the final state, conservatively assuming a common 
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a) pp -VW T + 2j, 




FIG. 6. Distribution of the invariant mass of the observable final state particles, m V i S , after 
all cuts, in pp —> + 2 j, £ / £', for the signal (solid line) with a) mg = 150 GeV and 

b) rriH = 180 GeV, and all backgrounds (except for the contributions from overlapping events 
and double parton scattering) at the VLHC (dashed: WWWjj, dotted: ttW, long-dashed: ttZ, 
long-dash-dot: ttj , dot-dashed: tttt). 
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FIG. 7. The m v i s distribution of the signal for pp —> £ ± £ /± + 4 j and mu = 180 GeV at a) 
the LHC, and b) the VLHC, in the SM (solid curve), for A hhh = A/A ski = 0 (dashed line) 
and for A hhh = 2 (dotted line). The dot-dashed line shows the combined m v i s distribution of 
all background processes (except those from overlapping events and double parton scattering). 
Qualitatively similar results are obtained for other values of mu- 
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lepton identification efficiency of e = 0.85 for each lepton. Except for the Higgs boson 
self-coupling we assume the SM to valid: by the time a measurement of A is attempted, 
the Higgs boson mass will be precisely known and the H —> W + W~ branching ratio will 
have been measured with a precision of 10% or better at the LHC and/or an e + e~ linear 
collider [45]. We include all background processes listed in Table II and Table III, except 
those from overlapping events and double parton scattering. The challenge of including 
higher order effects is considerably more complicated for the background than for the HH 
signal, where at least the physics interpretation is clear as previously discussed. The aim 
for the backgrounds is not to capture the bulk of events after cuts. Instead, one tries to 
cut into the tails of distributions, where the impact of higher order corrections might be 
very different. Therefore an analysis should depend as little as possible on the background 
rates [46], while a dependence on the signal rate is unfortunately unavoidable for any new 
physics process, which by definition will rely on comparably fewer, rare events. To show 
that our analysis fulfills this requirement, and approximately take into account the unknown 
NLO QCD effects, we perform two separate calculations of sensitivity limits: 

1. we assume a uniform K -factor of K = 1 for the rri ms distribution of the background 
but allow for a normalization uncertainty of AA f = 30% of the SM cross section; 


2. we assume a uniform A'-factor of K = 1.3 for the m V i S distribution of the background 
and allow for a normalization uncertainty of AA/" = 10% of the SM cross section. 


The results from both calculations are then compared and the more conservative bound is 
selected. Since the background cross section can be directly determined from the high rri ms 
region with a statistical precision of 15% or better for the assumed integrated luminosities, 
the bounds we derive should be conservative. 

The expression for \ 2 which we use to compute confidence levels is given by [47] 


X 


2 


E 


i= 1 


(Ni ~ fN ?) 2 
fN.r 


+ (- n D - 1 ), 


(17) 


where no is the number of bins, N t is the number of events for a given AA hhh = (A — 
A sm) / A sm, and A/ 0 is the number of events in the SM in the ith bin. / reflects the uncertainty 
in the normalization of the SM cross section within the allowed range, and is determined by 
minimizing y; 2 : 


with 


/ 


(1 + AA7)” 1 for / < (1 + AA/") _1 _ 

/ for (1 +AA7)" 1 < / < 1 +AA7 

1 + AA f for / > 1 + AA f 


P 



n D 

E 


i=l 



(18) 


(19) 


For the LHC, we derive sensitivity limits for integrated luminosities of 300 fb -1 , 600 fb -1 
and 3000 fb -1 , and Higgs boson masses in the range 150 GeV < < 200 GeV. An 

integrated luminosity of 300 fb -1 (600 fb^ 1 ) corresponds to 3 years of running at the LHC 
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FIG. 8. Limits achievable at 95% CL for AA hhh = (A — \sm)/Xsm in pp —> + 4 j at 

the LHC. Bounds are shown for integrated luminosities of 300 fb -1 (solid lines), 600 fb _1 (dashed 
lines) and 3000 fb -1 (dotted lines). The allowed region is between the two lines of equal texture. 
The Higgs boson self-coupling vanishes for A A hhh = — 1. 

design luminosity with one (two) detectors. The larger value of 3000 ftn 1 can be achieved in 
about 3 years of running at the SLHC with one detector. Since the cross section for Higgs 
boson pair production in the three lepton final state is very small, we calculate sensitivity 
bounds only for the same-sign dilepton channel. Our results are shown in Fig. 8, which 
demonstrates that, for 300 fb _1 , a vanishing Higgs boson self-coupling (AX HHH = —1) is 
excluded at the 95% CL or better, and that A can be determined with a precision of up to 
—60% and +200%. Doubling the integrated luminosity to 600 fb -1 improves the sensitivity 
by 10 — 25%. For 300 fb -1 and 600 fb -1 , the bounds for positive values of AA hhh are 
significantly weaker than those for A A hhh < 0, due to the limited number of signal events 
in this region of parameter space. At the SLHC, for 3000 fb -1 , the Higgs boson self-coupling 
can be determined with an accuracy of 20 — 30% for 160 GeV < rrin < 180 GeV. The 
significance of the SM signal for 300 fb _1 (3000 fb -1 ) is slightly more than 1 a (3 a) for 
rriff = 150 GeV and 200 GeV, and about 2.5 a (10 a) for Higgs boson masses between 
160 GeV and 180 GeV. The results shown in Fig. 8 are about 5 — 10% weaker than those 
found in Ref. [21] where only the dominant WWWjj and ttW backgrounds where taken into 
account while the effect of all other backgrounds was simulated by multiplying the combined 
WWW jj and ttW visible invariant mass distribution by a factor 1.1. 

For the VLHC, we calculate bounds for both the £ ± £ /± + 4 j and the f' ± C =f + ,=F + 2 j final 
states. We assume integrated luminosities of 300 fb -1 , 600 fb -1 and 1200 Do 1 . For a design 
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luminosity of C = 2 x 10 34 cm -2 s _1 [24], the latter corresponds to 3 years of running with 
two detectors. The 95% CL limits which one may hope to achieve at such a machine are 
shown in Fig. 9. 

At a pp collider with sfs = 200 TeV and an integrated luminosity of 300 fb _1 , the 
Higgs boson self-coupling can be measured with a precision of 8 — 25% at 95% CL for 
150 GeV < rriH < 200 GeV. For 1200 fb” 1 , the bounds improve to 4 — 11%. Although the 
signal to background ratio in the + 2 j channel is significantly better, the sensitivity 

limits which can be achieved are about a factor 2 to 3 weaker than those obtained for the 
+ 4 j final state, due to the reduced number of signal events. 

Our calculation of sensitivity bounds for A is subject to several uncertainties which should 
be addressed in a future more detailed analysis. In calculating limits, we have ignored 
the background from overlapping events and double parton scattering. Their sizes depend 
sensitively on the accelerator parameters, in particular the bunch spacing, and on the ability 
of the detectors to resolve the vertices of such events. In addition, these types of background 
are difficult to model at the parton level. Our estimates (see Sec. IIB 5) indicate that they 
are small at the LHC, but may not be negligible at the SLHC or VLHC. A conservative upper 
limit of how much the background from overlapping events and double parton scattering may 
change the bounds on A is obtained by assuming that the rn ms distribution of such events 
peaks in the region where also the signal reaches its maximum. Assuming that this is the 
case, and using the results of Table II, one finds that the SLHC (VLHC) limits weaken by 
at most 5% (15%). 

We also ignored the contributions from WWZjj and W + W~W + W~ production in our 
calculation. The cross section of these processes is small compared to that of the dominant 
background contributions. They therefore should have a negligible effect on the bounds 
which can be obtained. The extremely large number of Feynman diagrams contributing to 
WZjjjj and WWjjjj production makes a calculation employing exact matrix elements 
currently impractical. To calculate the cross section for these processes we interfaced the 
pp — > WZjj and pp — » WWjj matrix elements with PYTHIA. This procedure may well 
result in cross sections which differ from the correct result by a factor 1.5 ±:L . However, both 
the WZjjjj and the WWjjjj cross sections are small in the region of rn ms where the signal 
distribution peaks. Uncertainties in the calculation of their production cross sections thus 
should not change the bounds on A by more than a few per cent. To substantiate this claim, 
we have varied the WZjjjj cross section by a factor 2 ±1 and recomputed the 95% CL limits 
for A. The values obtained differ from those shown in Figs. 8 and 9 by at most 5%. 

Uncertainties in the extraction of sensitivity limits for A also arise from the ttj back¬ 
ground which we calculated at the parton level. Since the cross section of the ttj background 
sensitively depends on the lepton px cut and also the lepton isolation requirement, detector 
resolution effects may have a significant effect. Varying the ttj cross section by a factor 1.5 ±1 
changes the limits for A by about 5 — 10%. Finally, QCD corrections are expected to modify 
the shape of the m V i S distribution for both signal and background. In our calculation we 
have approximated the effect of QCD corrections by uniform AT-factors which do not take 
into account this effect. While an accurate answer how QCD corrections affect the shape 
of the m V i S distribution requires the calculation of next-to-leading corrections to signal and 
background processes, it seems unlikely that they will change the sensitivity bounds by more 
than 20%. 
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FIG. 9. Limits achievable at 95% CL for AA hhh = (A — Xsm)/Xsm for a) pp —■> + 4 j 

and b) pp —► £ ± £ /= F t' /=F + 2 j at the VLHC. Bounds are shown for integrated luminosities of 300 fb _1 
(solid lines), 600 fb -1 (dashed lines) and 1200 fb^ 1 (dotted lines). The allowed region is between 
the two lines of equal texture. The Higgs boson self-coupling vanishes for AA hhh = — 1 . 
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Our calculation of sensitivity bounds for A has been based on a simple y 2 test of the 
rn ms distribution. More powerful statistical tools, or a neural net analysis, may considerably 
improve the limits which can be achieved. 


V. DISCUSSION AND CONCLUSIONS 

A direct experimental investigation of the Higgs potential represents a conclusive test of 
the mechanism of clectroweak symmetry breaking and mass generation. After the discovery 
of an elementary Higgs boson and the test of its couplings to fermions and gauge bosons, 
experimental evidence that the shape of the Higgs potential has the form required for break¬ 
ing the clectroweak symmetry will complete the proof that the masses of fermions and weak 
bosons are generated by spontaneous symmetry breaking. In order to probe the shape of 
the Higgs potential, the Higgs boson self-coupling must be determined. 

The Higgs boson self-coupling can be measured in Higgs boson pair production at lepton 
or hadron colliders. In this paper, we presented a detailed analysis of Higgs boson pair 
production via gluon fusion with subsequent decay into four lU-bosons at the LHC, a lumi¬ 
nosity upgraded LHC (SLHC), and a planned next-generation hadron collider with a center 
of mass energy of y/s = 200 TeV (VLHC). We considered two final states: £ k l' ± + 4 jets and 
£ ± t' /=F £ /=F + 2 jets. To calculate the signal cross section, exact one-loop matrix elements for fi¬ 
nite top quark masses were used. Final state spin correlations for the H —> WW —> 4 fermion 
decay were fully taken into account, together with finite W and Higgs boson width effects. 

We investigated in detail which processes contribute to the background, including back¬ 
grounds from overlapping events and double parton scattering. All background cross sec¬ 
tions, except those for WWjjjj and WZjjjj production, were calculated using exact tree 
level matrix elements. Contributions to the background from overlapping events depend on 
the ability of detectors to resolve vertex positions, and on machine parameters. We pre¬ 
sented a simple order-of-magnitude estimate of the cross section from overlapping events 
which indicates that these should not be a problem at the LHC. At the SLHC and VLHC, 
however, the background from overlapping events could be non-negligible. 

At the LHC, the total background cross section is significantly larger than that of the 
signal in the £ ± £ /± -t-4j channel. There are too few events in the t' ± t' /=F t' /=F -|-2 j channel to make 
it useful. However, the distribution of the visible invariant mass of the final state particles, 
rn ms . for most of the processes contributing to the background peaks at a considerably higher 
value of rn ms than that of the signal, regardless of the value of A. The shape of the rn ms 
distribution can thus be used as a tool to derive limits on the Higgs boson self-coupling, A. 

At the VLHC, we found an improved signal to background ratio for the f' ± f" ± +4j channel. 
The t' ± t" =F t" =F + 2 j final state has an even more advantageous signal to background ratio, 
however, the signal cross section is significantly smaller than that for pp —> HH —> f' ± f" ± +4j. 

In order to determine how well one can hope to measure the Higgs boson self-coupling 
at future hadron colliders, we have performed a y 2 test of the rri ms distribution. We found 
that, at the LHC, with 300 flu 1 , one will be able to perform a first, albeit not very precise, 
measurement of the Higgs boson self-coupling. The non-vanishing of A, however, can be 
established at 95% CL or better for 150 GeV < rrin < 200 GeV. This alone is an important, 
non-trivial test of spontaneous symmetry breaking ; the exact non-zero value of A may vary 
depending on the way nature chooses to spontaneously break the electroweak symmetry. 
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At the SLHC, for 3000 fb -1 , a measurement with a precision of up to 20% at 95% CL is 
possible; A at the SLHC can be determined with an accuracy of 10 — 30% at the 1 a level 
for Higgs boson masses between 150 and 200 GeV. Compared with an estimate based on 
the total cross section [22], the fit to the rn ms distribution improves the accuracy of the 
measurement of Higgs self-coupling by a factor 1.2 to 2.5. For the same range of mn, the 
95% CL bounds on A for a 200 TeV pp collider (see Fig. 9) indicate that deviations of 10% 
or less from the SM value of A can be measured at 95% CL if more than 1 ab -1 can been 
accumulated. Due to the reduced signal rate, limits obtained from the £ ± F t £ /=f + 2 j final 
state are about a factor 2 to 3 weaker than those extracted from the + 4 j channel. 

It is interesting to compare the sensitivities which one may hope to achieve at the LHC, 
SLHC and VLHC with those obtained for future e + e~ linear colliders [15,17-19]. At TESLA 
energies, y/s = 500 — 800 GeV, the Higgs boson self-coupling can only be determined if 
rn h < 140 GeV. For larger values of mn, the cross section for the dominant Higgs pair 
production process, e + e _ —» ZHH , is too small for a useful measurement. For rn h = 
120 GeV, y/s = 500 GeV, and 1 ab -1 , one finds that A can be measured with a precision of 
5 A = ±0.20 for one sigma [18]. In contrast, Higgs boson pair production followed by decays 
into four W bosons at the LHC and SLHC offers an opportunity to probe the Higgs boson 
self-coupling for masses in the range 150 GeV < mn < 200 GeV. For mu < 140 GeV, where 
the decay H —> bb dominates, the QCD bbbb background is so large that a measurement of 
the Higgs boson self-coupling is hopeless. LHC and a linear collider operating in the range 
of y/s = 500 — 1000 GeV thus complement each other in their abilities to determine A. 

A more direct comparison can be carried out between the VLHC and CLIC, a proposed 
multi-TeV e + e _ linear collider [48]. For rn^ — 180 GeV, one finds [15] that, for e + e~ 
collisions at y/s = 3 TeV and an integrated luminosity of 5 ab -1 , A can be determined 
with a precision of 5 A = ±0.080 (la). For the same Higgs boson mass, the Higgs boson 
self-coupling can be measured with an accuracy of 5 A = ±0.035 at a 200 TeV pp collider 
with 300 fb -1 . 

Our analysis has been based on leading order parton level calculations. This introduces 
uncertainties in our derivation of sensitivity bounds which we estimated to be of 0(20%). 
In order to derive more realistic limits for the Higgs boson self-coupling, more detailed 
simulations which take into account detector effects, as well as the effects of higher order 
QCD corrections are needed. 
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